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The unbinding of vortex defects in the superconducting 
condensate with d-wave symmetry at T = is shown to lead 
to the insulator with incommensurate spin-density-wave or- 
der. The transition is similar to the spontaneous generation 
of the chiral mass in the three dimensional quantum electro- 
dynamics. Possible relation to recent experiments on under- 
doped cuprates is discussed. 

Common feature of all high-temperature superconduc- 
tors is that undoped they are Mott insulators with anti- 
ferromagnetic order [[[J. The central theme of the the- 
ories of cuprate superconductivity has therefore been 
to establish the connection between the insulating and 
the superconducting phase. Most of the work followed 
the usual route that suggests starting from the non- 
superconducting, in this case, Mott insulating phase, and 
trying to understand how it becomes superconducting. 
This approach was spectacularly successful for the con- 
ventional (low-T c ) superconductors, in part because the 
non-superconducting phase was a well understood metal- 
lic Fermi liquid. In cuprates, however, one does not enjoy 
this luxury, and the Mott insulator is strongly interact- 
ing and notoriously resistant to simple theoretical under- 
standing. This suggests one should look for alternative 
points of view that may be better adapted to the problem 
at hand. Since experimentally the superconducting phase 
seems to be a rather standard BCS-like d-wave state, one 
strategy would be to take this as a vantage point for fur- 
ther exploration of the cuprates phase diagram ||. 
Particularly interesting is the underdoped region, where 
experiments show a large pseudogap for spin excitations, 
and the superconductor-insulator transition at low tem- 
peratures. 

In this Letter I subscribe to the dual approach advo- 
cated above and show that the d-wave superconducting 
state (dSC) at T = has an instability towards the insu- 
lator with the incommensurate spin density wave (SDW) 
order. Using the Franz-Tesanovic transformation I de- 
rive the low-energy theory for the coupled system of d- 
wave quasiparticles and fluctuating vortices. Upon in- 
tegration over vortices the theory takes the form of the 
(anisotropic) 2+1 dimensional quantum electrodynamics 
(QED3) for two Dirac four-component spinors, which are 
related to the nodal quasiparticles by a singular gauge 
transformation, and are minimally coupled to the trans- 
verse gauge field ||. First, I show that the role of the 
coupling constant (or the "charge") in this gauge the- 
ory is at T =/= played by the thermodynamic fugacity 
of the vortex system. In the superconducting phase the 
"charge" is therefore zero, vortices are bound into pairs, 



and the gauge field is decoupled from the fermions. In 
the non-superconducting phase, on the other hand, the 
fugacity is finite, and the gauge field now mediates a long- 
range interaction between the Dirac fermions. The main 
result is that, at T — where the role of fugacity is played 
by the condensate of vortex loops, this interaction leads 
to an instability towards the incommensurate SDW or- 
der, through a condensed-matter equivalent of the chiral 
symmetry breaking phenomenon J4j. The T — tran- 
sition from the dSC into the SDW may be understood 
therefore as an instability of the gapless nodal fcrmionic 
excitation in presence of free topological defects towards 
the formation of bound states. Possible connections be- 
tween recent neutron scattering, ARPES, and STM ex- 
periments are discussed in light of this result. 

What follows rests on two postulates: 1) that there is a 
d-wave superconducting state in the phase diagram with 
sharp gapless quasiparticle excitations, and 2) that the 
amplitude of the superconducting order parameter may 
be assumed finite and inert much below the high pseudo- 
gap temperature T*, so that the only other relevant ex- 
citations are the topological defects in its phase (vortices 
and antivortices at T ^ 0, or vortex loops at T = 0). 
The first postulate is supported by the microwave ||, 
and the ARPES experiments ||, and the second by the 
measurements of the frequency dependent conductivity 
and the Nernst effect in the pseudogap regime ||. I 
begin by constructing the continuum, low-energy theory 
for the nodal quasiparticles in the d-wave state, using a 
different representation than in ]2j , QJ. The quasiparticle 
Hamiltonian is 

H qp = T ^2 i( iuj n~ £,j:)cl(k,uj n )c a (k,uj n ) (1) 

k,<T,UJ n 

-A(k)cl(k,u) n )cL a (-k, -w„) + c.c.], 

where A(fc) has the usual d-wave symmetry, and two 
spatial dimensions (2D) are assumed, c and are the 
electron operators, a = ± labels spin, and oj n are the 
fermionic Matsubara frequencies. In my units h = c = 
e = 1. Next, introduce two four- component Dirac spinors 

■^\ {2) {q,uj n ) = (c\{k,uj n ),C-(-k,-uj n ), (2) 

C+(fc - <3l(2),W„), C-(-k + Ql(2), -^n)), 

where Qin) = %Ki(2) ^ s the wavevector that connects 
the nodes within the diagonal pair 1(2). For the spinor 1, 
k = Ki+q, with \q\ <C \Kt\ (see Fig. 1), and analogously 
for the second pair. One has £V = and near the 

nodes. £r = — ff * . and Ar = -A,- * , for k f» 

' s fc S fe-Ql(2) ' k fe-Ql(2) ' 
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Ki(2)- Retaining only the low-energy modes in (1), and 
linearizing the spectrum as £g = Vfq x 
one arrives at the continuum field theory 



and = v A q y , 



S[9] = / d 2 r 



where = ^70 ■ 
02 <& 03 , and 72 = 



c!t*i[7o5 t + 7i«/<9 x + 72^^] *i + (3) 
(1 ->• 2,a;-> y,y-+ x), 

and the matrices 70 = 01 ® I, 71 = 
= — o~2 ® o\ satisfy the Clifford algebra 
Here a are the Pauli matrices, and the 



coordinate system has been rotated as in Fig. 1. 

Next, assume that the transition out of the dSC in 
the underdoped regime is due to unbinding of the topo- 
logical defects. This raises a rather non-trivial question 
of how to properly couple the vortex degrees of free- 
dom to quasiparticles Fortunately, this has recently 
been elegantly solved by Franz and Tesanovic 0, ||. 
Their idea is to split the phase of the order parame- 
ter A(r,r) = A| exp i(<p s (r, r) + <p r (f, r)), where s ( r ) 
is the singular (regular) part of the phase in presence of 
vortices, into two contributions, <p s + (f) r = 4>a + with 



N A 



V x V(/) A (r,T) = 2tt \ q lA S(r~ r iA {r)) 



iA=l 



(4) 



where qiA = ±1 is the unit vorticity of a (hc/2e) vor- 
tex (antivortex) defect, tracked by its coordinate in the 
imaginary time r*j(r), and analogously for vortices in B. 
Here V = (d T ,d x ,d y ). The division of vortices into two 
groups A and B is at this point arbitrary, and <p r is to 
be equally split between <j> A and 4>b- By making the 
singular gauge transformation in (1) from electrons into 
electrically neutral fermions c + r_\ — > c+(—) exp i4>A(B) > 
one immediately discovers that there is a hidden gauge 
field in the problem, a M = ^d^(4> A — 4>b), M = Oj 1) 2, that 
enters the theory (3) via minimal coupling <9 M — > d^ — ia^. 
The regular part of the phase 4> r cancels in a, which is 
entirely due to vortices. <j) r is contained in the second, 
Doppler gauge field, t> M = \d^(4> A + 4>b), which enters 
the theory for the neutral fermions precisely as the true 
electromagnetic gauge field would. Gauge invariance pro- 
tects a from becoming massive from the integration over 
fermions, while such a protectorate does not exist for v. 
Power counting implies then that the coupling of fermions 
to v is irrelevant, and may and will therefore be dropped 
hereafter. 

It was argued in || that although a can not become 
massive from fermions, it should be massive if vortices are 
bound. Next I present a simple but a rigorous derivation 
of the dynamics of the gauge-field a at T ^ 0, where one 
can get away with the neglect of the quantum fluctua- 
tions, which supports this insight. Assume a collection 
of A + (iV_) vortices (antivortices) at the positions {u}. 
The energy of the (classical) vortex system is 



1 N 



rj), 



(5) 
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FIG. 1. The wavevectors Ki, K2, and q. 



where v(r) ks — In \r\, at large distances, and N = N + + 
N~. The partition function of the coupled system of 
quasiparticles and vortices can be written as 



Z = J D[^]e~ S[<S! ^Z v , 



(6) 



where Z v is the grand-canonical classical partition func- 
tion of the vortex system (2D Coulomb plasma) 



Z v — 



E 
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N-,N+N-,N+ A B A B 



(7) 



where = iVt ( ) + 



Ng ', and y is the bare vor- 



tex fugacity. To preserve the a — > —a symmetry in the 
original Hamiltonian in Z v I average over all possible di- 
visions of vortices and antivortices into groups A and B. 
This ensures that on average there is an equal number of 
vortices (and antivortices) in both groups. Next, intro- 
duce the vorticity densities in Z v by inserting the unity 



r N A 

1=1 D[p A ]5(p A (r) - ^ qiA$(r- r iA )) 



(8) 



and similarly for B. The gauge field then becomes 

(V x a{r)) T = 2ir(p A (r) - p B (f)), (9) 

in the transverse gauge V • a = 0. Subindex r de- 
notes the r-component. v is defined the same way ex- 
cept with the plus sign between p A and ps- Perform- 
ing then the Gaussian integrations over p A , ps, and 
v, and the summations in the Eq. (7) exactly, yields 
Z v = J £>[$+,$_, a] exp(-S„[$ + ,$_,a]), with 
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S v [$ + ,$-,a\ = J d 2 f\2T(V<P + (r)) 2 + 
$_(f)(V x d(r)) T - 2ycos(3> + (r)) cos($_(r))]. 



(10) 



Real fields $± = $ A ± $b are the Lagrange multipliers 
introduced to enforce the constraints in the Eq. (8) p0|. 
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The partition function of the coupled system of d-wave 
quasiparticles and vortices at T ^ is therefore Z = 
/£>[*,$+,$_, a] exp(-5[*, a] - S v [$+,$_, a]) with 

F r r 13 

S[V, S] = V / d 2 f / drf - (11) 

<=i J J ° 

with F = 2, and the x < > y exchange of the coordinates 

for the i = 2 component is assumed. I have also set Vf = 
i>a = 1 here for simplicity. The Dirac field \& represents 
the neutral (gauge-transformed) fermions, and S v is given 
by the Eq. (10). This is my first result. It has several 
remarkable features. First, if one turns off the coupling 
to fermions (by taking, formally, the quenched limit F = 
0), the integration over a in (10) simply enforces = 
0. S v reduces then to the standard sine-Gordon theory, 
which is known to provide the correct description of the 
Kosterlitz-Thouless transition |ll[] . More importantly for 
our purposes, for F = one also finds 



((V x a(f)) T (V x 3(f)) T ) = (y)S(r- j*), 



(12) 



where (y) = y(27r) 2 (exp(i$ + )}, with the average to be 
taken over S v with $_ = 0. (y) may be recognized as 
the thermodynamic, or the renormalized, fugacity of the 
vortex system [Q. The integration over the fields $+ 
and $_ is thus equivalent to reducing the action (10) to 



S v 



2(y> 



(13) 



in the partition function Z . In the dielectric phase of the 
vortex system the field <£>+ is massless, and consequently 
(y) = fPH , so the gauge- field asymptotically decou- 
ples from the fermions. Quasi-particles become sharp 
excitations in the dSC, in agreement with the ARPES 
H and the microwave measurements [^). In the non- 
superconducting phase, on the other hand, vortices are 
free, becomes massive, and (y) ^ 0. This has pro- 
found consequences for the fermions, as I discuss shortly. 

At T — quantum fluctuations need to be included, 
as the topological defects in 2 + 1 dimensions become 
vortex loops |12[| . It seems clear on physical grounds, 
however, that after the integration over the loops, (apart 
from the inherent anisotropy,) the form of the action for 
the gauge field should remain similar to the Eq. (13), 
except that (V x a) 2 T /(2{y)) -> (V x a) 2 /{2(y)). This 
can also be derived on a lattice, where one finds that the 
role of the coupling (y) at T = is assumed by the dual 
order parameter that becomes finite only when there are 
infinitely large loops in the system, and which is tanta- 
mount to the loss of phase coherence . This way one 
finally arrives at the QED3 with the full Maxwell term 
for the transverse gauge-field a as the relevant low-energy 
theory. 

QED3 has been extensively studied by field theorists 
as a non-trivial toy model exhibiting the phenomena of 
dynamical symmetry breaking and confinement B. In 



particular, it has been established that for the number of 
Dirac fields F < F c the interaction with the gauge-field 
is strong enough to spontaneously generate the so-called 
chiral mass for fermions. I will demonstrate that the chi- 
ral mass in the theory (11) is nothing but the SDW order 
parameter. First, to acquire some sense for the chiral in- 
stability consider the fermion propagator. Neglecting the 
vertex and the wave-function renormalizations, it can be 
written as G~ 1 (p) = i^f v p v + S(p), where the self-energy 
satisfies the self-consistent equation 



(2tt) 3 p 2 + Y?(p) 



(14) 



with q = (Lo,q x ,q y ). The gauge-field propagator in the 
transverse gauge is — {8^ v ~p^p v )/(p 2 +U(p)), where 
II(p) is the self-consistently determined polarization. As- 
suming m = S(0) 7^ gives Q 



n(g) 



(y)F 

2ir 



(m 



Am 2 



2q 



\/ q 2 + 4m 2 ' 



(15) 



When this is inserted into the Eq. (14), it can be shown 
that there is a solution with a finite m only when F < F c , 
with F c — 32 /it 2 Q. More elaborate calculations that 
fully include the wave-function renormalization and the 
vertex corrections confirm this result, and yield F c k, 3 
pTj . Simulations on the lattice version of the QED3 g| 
also find 3 < F c < 4, in agreement with the analytical 
estimates. 

By reversing the transformations that led to the QED3 
the reader can convince himself that the mass term for 
neutral fermions is equivalent to the low-energy part of 
the following addition to the electronic Hamiltonian (1): 



mT 



oc\{k + q,u> n )c a (k,uj n ), 



(16) 



so that the chiral mass may be identified with the SDW 
order parameter (or the staggered potential) along the 
spin z-axis, and at the wavevectors Qi,2- This is, of 
course, why the particular construction of the Dirac field 
was made in the first place. Note that: 1) the SDW or- 
der is induced already at an infinitesimal vortex fugacity, 
but it is rather weak, rn « (y) / exp(27r / yj (F c / F) — 1) 
for F ks F c Q, and 2) neutral fermions are bound (con- 
fined) at large distances in the SDW, by the weak log- 
arithmic potential (provided by the fact that 11(g) = 
(y)Fq 2 j '((mm) for q <C rn). With some anisotropy 
(vf ^ va), the global symmetry of the massless theory is 
only U(2) x U(2), so the mass term reduces each U(2) to 
U(l) x U(l). The two broken generators per Dirac field 
rotate the "cos-SDW" in the Eq. (16) into either the sim- 
ilar "sin-SDW", or into the phase- incoherent state with 
"d + ip v pairing between the neutral fermions. SDW is 
therefore the only state that can be obtained by the un- 
binding of vortex defects and that respects parity. In 
the isotropic limit the massless theory recovers the full 
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U(4) symmetry and additional broken symmetry states 
become available, like, remarkably, the stripe-like charge 
density waves parallel to a or b axis [[l3| . 

Once it is realized that unbinding of vortices leads to 
the SDW order, it becomes natural to wonder what the 
nature of vortices inside the dSC could be. From the per- 
spective of this work it seems more than plausible that 
vortex cores are actually in the insulating phase, so that 
by approaching half-filling one lowers the core energy. In 
this picture the chemical potential should be related to 
the bare vortex fugacity, which when too large leads to 
the proliferation of defects, in analogy to the Berezinskii- 
Kosterlitz-Thouless transition |ll|] . The idea of SDW in 
vortex cores finds some experimental support in the re- 
cent STM [H| and the neutron scattering studies [^7| , 
as well as in the mean-field calculations Eq). Further- 
more, the present work suggests that the superconductor- 
insulator transition should be accompanied by the ap- 
pearance of the incommensurate SDW correlations at 
the wave vectors Qi,2, with the incommensurability in- 
creasing with doping. This is consistent with the recent 
neutron scattering experiments |l9[ ] on the underdoped 
LaSrCuO close to the superconducting transition. Fi- 
nally, the d-wave pseudogap should continuously evolve 
into the insulating state, except for the gap that should 
develop at the nodes. This also seems in agreement with 
the observations po[ . 

To summarize, I showed how liberating topological de- 
fects in the d-wave superconductor at T = leads to the 
incommensurate SDW, which is then expected to con- 
tinuously evolve into the commensurate antifcrromagnct 
close to half-filling. Near the transition the SDW order 
is inherently weak due to the relative closeness of the two 
flavor QED3 to its chiral critical point at F c w 3. The 
SDW transition temperature near the superconductor- 
insulator transition may therefore be expected to be 
much lower than the corresponding superconducting T c 
on the other side of the transition, not in contradiction 
with the known topology of the cuprates phase diagram. 
The issue of a quantum-disordered (deconfined) ground 
state in this approach reduces to whether F c may in fact 
be smaller than two. It has recently been argued that 
F c < 3/2 in spite of the results of virtually ev- 

ery calculation violating this bound. Recent numerical 
study on larger systems than in fl5|j , for example, find a 
small, but definitely a finite mass at F = 2 p2fl . More 
interesting possibility is that F c may depend on some ad- 
ditional parameter in the theory. For example, one may 
speculate that a large anisotropy Vf/vA ~ 10, which is 
certainly present in cuprates, could affect F c . This way 
F c would become doping dependent, which could open a 
route for the T = deconfined phase in between the dSC 
and the SDW. Further examination of this possibility, 
together with the issue of the order of transition, effects 
of disorder and finite temperature, and the relations to 
other theoretical approaches to the high-T c problem will 
be discussed in the future longer publication |13[ . 
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